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Abstract. We present a detailed phenomenological study of spin-dependent single-inclusive high-pr hadron 
photoproduction with particular emphasis on the kinematics relevant for the COMPASS and Hermes fixed- 
target experiments. We carefully examine the theoretical uncertainties associated with the only moderate 
transverse momenta accessible in such measurements and analyze the sensitivity of the relevant spin 
asymmetries to the gluon polarization in the nucleon as well as to the completely unknown parton content 
of circularly polarized photons. 

PACS. 13.88. +e - 12.38.Bx - 13.85.Ni 



1 Motivation and Introduction 

To measure the so far largely unknown polarization of 
gluons, Ag, in the nucleon is a key goal of several cur- 
rent experiments in high-energy nuclear physics. The suc- 
cessful start of polarized proton-proton collisions at the 
BNL Relativistic Heavy Ion Collider (RHIC) has opened 
up new, unequaled possibilities 1 . Gluon polarization can 
be accessed in a variety of high transverse momentum pro- 
cesses such as single- inclusive hadron, jet, prompt photon, 
or heavy quark production. Compass j2J at CERN and 
Hermes |3j at DESY, instead, scatter beams of longitu- 
dinally polarized leptons I off polarized fixed targets N. 
Here, in particular high-py hadron pairs, both in photo- 
production and in deep-inelastic electroproduction, have 
been identified by the experiments to provide a promising 
tool to gain some knowledge about Ag 0|. 

A meaningful extraction of parton densities from ex- 
periment requires a reliable interpretation of the underly- 
ing data in the framework of perturbative QCD (pQCD). 
Factorization theorems [5] assure that in the presence of 
a large scale, for instance, a high transverse momentum 
Pt, a cross section can be written as a convolution of par- 
tonic hard scattering cross sections, which depend on the 
process under consideration, certain combinations of non- 
perturbative, universal parton densities, and, if applicable, 
fragmentation functions. The standard pQCD framework 
receives corrections which are suppressed by inverse pow- 
ers of the large scale characterizing the process. Since it 
is not a priori clear where these "power corrections" be- 
come relevant, the factorized ansatz has to be scrutinized 
for each measurement. For all spin-dependent observables, 
cross sections or spin asymmetries, this is best achieved 



by studying first the unpolarized cross sections because 
both ingredients, the parton densities and the hard scat- 
tering cross sections, are known very well here. Only if a 
reasonable agreement between theory and experiment is 
established within their respective uncertainties can one 
have confidence that a similar measurement with polar- 
ization can be used to extract information about the spin 
structure of the nucleon. 

First measurements at RHIC of wnpolarized high-pT 
neutral pion [S] and prompt photon [7] production rates 
turned out to be in good agreement with theoretical pre- 
dictions 8,9, : TO] even down to unexpectedly small values 
of the transverse momentum of about 1.5 GeV. For this 
comparison it is essential to compute the perturbatively 
calculable short-distance cross sections up to the next- 
to-leading order (NLO) approximation of pQCD. On the 
one hand, NLO corrections turn out to be sizable, on the 
other hand, theoretical uncertainties related to the renor- 
malization and factorization procedures are much reduced 
compared to lowest order (LO) estimates. Therefore much 
effort was put in calculating the NLO pQCD corrections 
also for spin-dependent processes in recent years. 

The Phenix collaboration at RHIC has recently re- 
leased first, partly preliminary, results for high-p^ neu- 
tral pion production with longitudinally polarized protons 
|TTj . However, from this measurement alone it is difficult 
and way too early to draw any conclusions about the gluon 
polarization Ag in the nucleon; see also [12] . To map out 
Ag(x,/j.) over a sufficiently large range in the momentum 
fraction x, which is ultimately required for a determina- 
tion of the first moment of Ag, and scale /it (to test the 
scale evolution predicted by pQCD) a much better statis- 
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tical accuracy and a significantly larger coverage in px are 
needed. As a rule of thumb, a measurement of a final state 
jet or hadron with a certain px predominantly probes the 
gluon density at momentum fractions x > 2px/\^S where 
y/S is the available proton-proton center-of-mass system 
(c.m.s.) energy. Since the cross section drops rapidly with 
increasing pt it will take RHIC quite some time to accu- 
mulate enough events to determine Ag at x > 0.1. This 
is where the fixed-target experiments at much lower en- 
ergies can possibly help to add information about Ag at 
medium-to-large x in the very near future. 

In the case of hadron-pair production, IN — ► VH1H2X, 
the process already considered by Compass and Hermes 
the theoretical framework is rather complex and so 
far NLO corrections are still lacking. In this paper we 
therefore focus on the photoproduction of single-inclusive 
hadrons, IN — > I'HX, where the complete NLO pQCD 
framework is already at hand [51ll3| and has been recently 
applied to study the physics case of a conceivable first 
polarized ep-collider at BNL j!4j . Photoproduction with 
quasi-real photons has the advantage of much higher rates 
than deeply-inelastic electroproduction of hadrons. The 
price to pay is the more involved theoretical framework for 
photoproduction where so-called "direct" and "resolved" 
photons contribute to the cross section as depicted in 
Figs. E (a) and (b), respectively. Estimates for the latter 
require knowledge of the parton content of circularly po- 
larized photons which is lacking at the moment. We will 
demonstrate below that this does not, however, seriously 
limit the usefulness of this process. 

We note that there are actually already some photo- 
production data for double-longitudinal spin asymmetries 
in eiV -> H±X from the El 55 experiment at SLAC 
Here the scattered electron was not observed, so the data 
set is integrated over all virtualities of the exchanged pho- 
ton. Nonetheless, thanks to the photon propagator, pro- 
duction by almost real photons strongly dominates. The 
El 55 data have been compared to the theoretical predic- 
tions of 53 , which were based on LO calculations supple- 
mented by certain higher-twist contributions. A relatively 
poor overall agreement was found. To achieve a better 
description of the data, the theoretical framework used 
in J7| was subsequently augmented by estimates JS| of 
soft contributions that are mostly relevant at small trans- 



verse momenta. We will later present a brief comparison 
of our NLO calculations with the El 55 data. 

At fixed-target energies it is much less clear that pQCD 
methods are applicable than at collider energies. In fact, 
data for high-py processes in hadron-hadron fixed-target 
scattering have in the past been a serious challenge for 
the factorized framework outlined above ^Hj- It is there- 
fore even more important to demonstrate first that unpo- 
larized production rates can be described with standard 
pQCD methods before turning to studies of spin asym- 
metries. Otherwise conclusions about the parton and/or 
spin content of nucleons might be misleading or wrong. We 
will therefore not only present results for spin asymmetries 
and discuss their possible sensitivity to Ag but also focus 
on predictions for unpolarized reference or "benchmark" 
cross sections which allow to probe the pQCD framework. 

The paper is organized as follows: in Sec.|2]we briefly 
recall the theoretical framework for photoproduction re- 
quired for the subsequent numerical sections and set up 
our notation. We refrain from repeating lengthy technical 
details and refer the interested reader to Refs. |%1TH5] . We 
then turn to a detailed study of single-inclusive pion pro- 
duction in the kinematical regions relevant for COMPASS 
and Hermes including, as far as possible, experimental 
cuts. Detailed phenomenological results are presented in 
Sees. 13. II and 13. 21 for COMPASS and Hermes, respectively. 
Special emphasis is put on theoretical uncertainties asso- 
ciated with the application of perturbative methods in the 
domain of the rather moderate transverse momenta acces- 
sible in fixed-target experiments. We first present unpolar- 
ized "benchmark" cross sections for each experiment and 
then turn to a discussion of the relevant spin asymmetries 
and their sensitivity to Ag. We comment on the relevance 
of resolved photon contributions and their influence on the 
extraction of the gluon polarization. In Sec. 13. 31 we present 
the comparison to the E155 data. We conclude in Sec.0] 



2 Technical Prerequisites 

In the following, we consider the spin-dependent photo- 
production cross section IN — > I'HX for the scattering of 
a longitudinally polarized lepton beam I off a target nu- 
cleon N with subsequent hadronization of a produced par- 
ton c into the observed hadron H with momentum Ph, see 
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Fig. ^ The observed hadron H is at high transverse mo- 
mentum px, ensuring the required large momentum trans- 
fer. The differential single-inclusive cross section can then 
be schematically written as a convolution 



dAa — — [da ++ — da^ 



(1) 



= ^ / dx a dx b dz c Af l a (x a ,[i f )Af b N (x b ,[if) 

,-, I, „ J 



xdAa ab ^ cX (S, x a , x b , P H / z c , [i r , [i f , [i' f )D c {z c ,fi' f ) . 

(2) 

The subscripts "++" and "H — " in Eq. denote the 
helicities of the colliding leptons and nucleons. 

In Eq. J5J x b is the fraction of the nucleon momentum 
taken by parton b, and the Af^ (x b , [if) are the usual po- 
larized parton densities. For instance, the spin-dependent 
gluon distribution at a scale [if is defined by 



Ag (x b , [if) = g+(x b , [if) ~ g_ (x b , [if) 



(3) 



where the subscripts ± indicate the helicity of a gluon in a 
nucleon N of positive helicity. The non-perturbative func- 
tions Df? (z c , n'j) describe the fragmentation of parton c 
(c = q, q, g) into the observed hadron H at a momentum 
scale [if. z c is the fraction of parton c's momentum taken 
by the hadron H . The sum in Eq. (0) runs over all par- 
tonic channels a + b — > c + X contributing to the single- 
inclusive cross section IN — > I'HX with dAa ab -* c x the 
associated spin-dependent partonic hard scattering cross 
sections. The latter can be calculated in pQCD order-by- 
order in the strong coupling constant a s (ii r ), with [i r de- 
noting the renormalization scale, and are known up to 
NLO accuracy. The spin-dependent results can be found 
in [51ll3U2(J| . The experimentally measured cross section is 
the sum of the so-called "direct" and "resolved" photon 
contributions, cf. Figs. Q (a) an d (b), respectively, 



dAa = dAa^ir + dAa rc 



(4) 



which both can be cast into the form of Eqs. (JTJ and J5J 
by defining Af l a appropriately. Most generally, the parton 
densities of the lepton I, Af l a: can be written as convolu- 
tions, 



upper limit Q^ax (determined by the experimental condi- 
tions) by a lepton of mass to/. 

The explicit form of the Af2 in Eq. (JjjJ depends on 
the specifics of the interaction that the quasi-real photon 
radiated off the lepton undergoes in the hard scattering. 
For the direct photon contribution to the cross section, 
dAa^h-, depicted in Fig. ^ (a), parton a in Eq. © has 
to be identified with an elementary photon and hence x a 
with the momentum fraction y of the photon w.r.t. the 
parent lepton, i.e., 



Af2 = 5(1 - x,) 



(7) 



in Eq. © . If the photon first resolves into hadronic "con- 
stituents" which then interact with partons of the target 
nucleon N, Af2 in Eq. JSJ represents the parton densi- 
ties of a circularly polarized photon. The latter are de- 
fined in complete analogy to the ones for a nucleon tar- 
get in Eq. Q . Unlike hadronic parton distributions which 
are genuinely non-perturbative objects, photonic densities 
consist of a perturbatively calculable "pointlike" contribu- 
tion which dominates their behavior at large momentum 
fractions x 1 in Eq. (jSJ . This will become important in the 
discussion of the numerical results in the remainder of the 
paper. 

We shall note that neither dAa^ nor dAa rcs indi- 
vidually are measurable cross sections as their definition 
depends on the specific details of the factorization proce- 
dure for dealing with singular collinear parton emissions. 
Only their sum, Eq. (@J, is a physically meaningful quan- 
tity provided both dAa^ir and dAa rcs have been calcu- 
lated consistently in the same factorization scheme. Need- 
less to say, the corresponding spin-averaged cross sections 
are straightforwardly obtained by replacing all polarized 
parton densities, partonic cross sections, etc., by their ap- 
propriate unpolarized counterparts. 

Finally, we wish to stress that Eqs. (JTJ and J5J apply 
to a single-inclusive cross section, which is not to be con- 
fused with a cross section for producing a so-called "lead- 
ing hadron". The latter, for which only the hadron with 
the highest momentum in the event is counted, is some- 
times provided by experiment. It is not possible in the 
framework of pQCD based on the factorization theorem 
and outlined above to compute a "leading-hadron" cross 
section, as fragmentation functions Df? are fully inclusive. 



Af a (x a>t if) = f ^ AP^y) Af2 (x 

j x a y v 



with 



i - (i - vf 



hi- 



<3L(! - v) 



2 9 

m iV 



+ Zrnfy 2 



) 2 

J max 



i-t/ s 

mfy 2 / 



(5) 



(6) 



being the spin-dependent Weizsacker- Williams "equiva- 
lent photon" spectrum |2J that describes the collinear 
emission of a photon with low virtuality Q 2 less than some 



3 Phenomenological Results 

With the technical framework and notation set up we are 
now in a position to perform a detailed phenomenologi- 
cal analysis of polarized and unpolarized single-inclusive 
pion photoproduction cross sections for the Compass and 
Hermes fixed-target experiments in Sees. 13. ll and 13.21 re- 
spectively. In Sec. 3.3 we present the comparison to the 
E155 data. 

All our results will be differential in the hadron's trans- 
verse momentum px and integrated over all kinematically 
and experimentally allowed pseudo-rapidities t] of the pro- 
duced hadron H measured w.r.t. the direction of the inci- 
dent lepton beam. For the unpolarized parton densities of 



4 



B. Jager et al.: Longitudinally Polarized Photoproduction of Inclusive Hadrons at Fixed- Target Experiments 



the nucleoli and photon we adopt throughout the CTEQ 
[22] and GRV (2H| sets, respectively. To study the sensi- 
tivity to the unknown gluon polarization of the nucleon 
we use four different sets of spin-dependent parton dis- 
tributions emerging from the GRSV analysis [21] • These 
sets span a rather large range of gluon densities Ag all 
consistent with present DIS data. Apart from the "stan- 
dard" set of GRSV with a moderately large, positive Ag, 
the three other sets "Ag = g input", "Ag = input", 
and "Ag = —g input" are characterized by a large posi- 
tive, a vanishing, and a large negative gluon polarization, 
respectively, at the input scale of the evolution. 

To model the hadronization into pions, we use the frag- 
mentation functions of KKP |25| which are known to pro- 
vide a good description of the pp — > ir°X cross sections 
measured by Phenix and Star [S] (and, of course, all 
available e + e~ data). Although the framework outlined 
in Sec. [21 applies in principle to the photoproduction of 
any hadron species (made of light quark flavors only) , we 
limit ourselves mainly to the production of neutral pions. 
These are most abundant, and the tt° fragmentation func- 
tions, mainly extracted at large resolution scales fi ~ Mz 
from LEP data, have been shown to work quite well also 
at scales of a few GeV relevant here |26| . 

Since the parton distributions of circularly polarized 
photons, Ap , required for estimates of the resolved pho- 
ton contribution <iZ\cr rcs , are completely unknown so far 
one has to resort to some models. For our purposes it is 
sufficient to use the two extreme scenarios proposed in 
Ref. [2H| which are based on "maximal" [Af(x, /Lto) = 
P(x,/j,q)] and "minimal" [Ap(x, /xq) = 0] saturation of 
the positivity condition \Ap 1 (x, fio)\ < p(x,[io) at the 
starting scale /iq for the evolution to /i > /io . Both models 
result in very different parton distributions Ap at small- 
to-medium £ 7 while they almost coincide as x 7 — > 1 due to 
the dominance of the perturbatively calculable pointlikc 
contribution in this region. Unless we study the impact of 
the details of the non-perturbative input to the evolution 
of Ap on the the full photoproduction cross section, the 
use of the "maximal" set will be implicitly understood. 
For recent work on the parton distributions of polarized 
photons, see |27j. 



3.1 Single-Inclusive Pion Production at COMPASS 

With the present setup, Compass scatters polarized muons 
off the deuterons in a polarized 6 LiD solid-state target. 
The beam energy is = 160 GeV, corresponding to a 
lepton-nucleon c.m.s. energy of ^/~S ~ 18 GeV. On aver- 
age the beam polarization is ~ 76%. About — 50% 
( "dilution factor" ) of the deuterons can be polarized, with 
an average polarization of Vd — 50% [2] ■ It is also conceiv- 
able to run in the future with a "proton target" (NH3) 
which, however, will have a much less advantageous di- 
lution factor of about T p ~ 17.6% but a polarization of 
V p ~ 85%. We will therefore mainly present results for the 
photoproduction of neutral pions off a deuteron target. 

Currently, pions can be detected if their scattering an- 
gle is less than 6> max = 70 mrad. Using 77 = — In tan (0/2), 
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Fig. 2. Unpolarized and polarized pT-differential single- 
inclusive cross sections at LO (dashed) and NLO (solid) for 
the photoproduction of neutral pions, pud —* fj,'n "X at jS = 
18 GeV, integrated over the angular acceptance of COMPASS. 
The lower panel shows the ratios of NLO to LO contributions 
(A"- factor). 



this straightforwardly translates into a minimal bound on 
pseudo-rapidity n nlln ~ 3.35 in the laboratory frame, cor- 
responding to — 0.44 in the lepton-nucleon c.m.s. 
where we have made use of the well-known behavior of 
rapidity under Lorentz boosts: 



1. 



(8) 



with Mm the nucleon mass. Kinematics sets an additional 
upper bound on the pion's rapidity depending on its trans- 
verse momentum: ryJ^J = cosh _1 (v / 5'/2pT). As already 
mentioned, we will always integrate over all kinematically 
allowed rapidities, 0.44 < rj cms < cosh _1 (\/>S'/2pT)- 

The momentum distribution of the quasi-real photons 
radiated off the muons can be described by the Weiz- 
sacker- Williams spectrum given in Eq. © with mi — m M 
and Q^ ax = 0.5 GeV 2 . The photon's momentum frac- 
tion y is restricted to be in the range 0.2 < y < 0.9. At 
smaller y the photon polarization is strongly diluted as 
the unpolarized equivalent photon spectrum behaves like 
(1 + (1 - yf)/y rather than (1 - (1 - y) 2 )/y in Eq. ©. 
We note that the often omitted non-logarithmic pieces in 
Eq. jnj result in a small but non-negligible contribution 
here. 

Figure [21 shows our predictions for the p^-differential 
cross section for the polarized and unpolarized photopro- 
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Fig. 3. Scale dependence of the LO and NLO unpolarized 
(a) and polarized (b) cross sections for fid — > fi'i^X shown 
in Fig. |21 All scales are varied simultaneously in the range 
Pt/2 < fi r = fj,f = n'f < 2pr- Solid lines correspond to the 
choice where all scales are set to pr- All LO computations 
have been rescaled by a factor 0.1 to better distinguish them 
from the NLO results. 



duction of neutral pions, fid — > fi'ir°X, at LO and NLO 
accuracy at Compass. We have set all factorization and 
renormalization scales in Eq. equal to the pion trans- 
verse momentum, fif = fi'f = fx r = pt- The so-called 
i\-factor, defined as the ratio of NLO to LO unpolarized 
(polarized) cross sections, 



K 



d(A)a 



NLO 



d{A)a 



LO 



(9) 



is depicted in the lower panel of Fig. [5J The results indi- 
cate clearly the relevance of the NLO corrections to the 
cross section in the small-to-medium p^-region accessible 
at a fixed-target experiment, in particular in the unpo- 
larized case. The effects of higher order corrections are 
less pronounced for the longitudinally polarized cross sec- 
tion, mainly due to large cancellations among the various 
partonic channels contributing to dAa at NLO. The dif- 
ferent behavior of the unpolarized and polarized if-factors 
also clearly indicates that the contributions due to NLO 
corrections do not cancel in the experimentally relevant 
double-spin asymmetry defined as 



.4 



H,N _ 
LL — 



dAu du-L 



dcr -i 



da (icTa 



da j 



(10) 



Not unexpectedly, at pr — 1^2 GeV, the K- factors for 
both, polarized and unpolarized cross sections rise sharply, 
perhaps indicating a breakdown of the standard pQCD 
framework as outlined in Sec. [21 Since, as will be demon- 
strated below, this is precisely the p-r-region where the 
statistical accuracy of Compass would best allow to de- 
duce some information about Ag from a measurement of 
the double-spin asymmetry A^jf, one has to ensure the 
validity of the pQCD framework first. As already empha- 
sized in the Introduction, this is best achieved by a mea- 
surement of the unpolarized cross section shown in Fig. [2 
where all ingredients, partonic cross sections, parton dis- 
tributions, and fragmentation functions are known. We 
note that all-order resummations of large logarithms in 
the perturbative series which appear when the initial par- 
tons have just enough energy to produce a high-px pion 
and a recoiling massless "jet" may lead to a considerable 
enhancement of the cross section at fixed-target energies 
as was recently demonstrated for the process pp — > ttX 
|29j . Similar calculations for the case of photoproduction 
are not yet available but certainly desirable. Any resid- 
ual shortfall of the rcsummed theoretical prediction would 
then indicate the relevance of non-perturbative contribu- 
tions. 

Large X-factors, as found in Fig. [3 are, however, of 
limited significance for unambiguously estimating the im- 
pact of higher order corrections in a perturbative calcula- 
tion. This is due to the large scale uncertainties associated 
with the LO cross sections entering the denominator of 
Eq. We therefore further explore the reliability of the 
perturbative approach by studying the dependence of the 
calculated cross sections, Eq. on the unphysical, a pri- 
ori arbitrary factorization and renormalization scales, fif, 
fi'j and fi r , respectively. Any dependence on these scales 
is a remnant of the truncation of the perturbation series 
at some fixed order of a s and thus expected to dimin- 
ish if higher order corrections are included. This is the 
prime motivation for going beyond the LO approximation 
of pQCD. The scales are of the order of the hard scale char- 
acterizing the process, here, the large px of the observed 
hadron, but not further specified by theory. An estimate 
for the sensitivity of the computed cross section to fif, p/p 
and fi r is usually obtained by varying them collectively in 
the range pt/2 < fif = fi'f = fi r < 2pr- We note that in 
principle all scales can be varied independently. 

The shaded bands in Figs. 01 (a) and (b) indicate the re- 
sulting scale uncertainty of the unpolarized and polarized 
cross sections, respectively, shown in Fig.|3 In contrast to 
similar studies for inclusive high-py pion hadro- [B] and 
photoproduction 13 at collider energies, where the theo- 
retical scale uncertainties are substantially reduced when 
the NLO corrections are taken into account, this barely 
happens here. This is particularly true for the unpolar- 
ized cross section, whereas the scale dependence of the 
polarized cross section improves beyond the LO, but only 
slightly in the region 1 < px < 2 GeV which mainly mat- 
ters for a determination of Ag. Together with the large 
if-factors found in Fig. |21 this underlines the delicacy of 
a perturbative calculation in the low-energy range associ- 
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Fig. 4. Double-spin asymmetry at NLO for different 

gluon polarizations in the nucleon and minimal and maximal 
saturation of the polarized photon densities, dashed and solid 
lines, respectively (see text). The "error bars" indicate the esti- 
mated statistical uncertainty for such a measurement at COM- 
PASS. 



ated with fixed-target experiments such as COMPASS. It is 
therefore particularly important to check the applicabil- 
ity of pQCD methods by showing, for instance, that data 
taken in unpolarized collisions fall within the uncertainty 
band shown in Fig. 

Next we consider the double spin asymmetry A^jf, 
Eq. I|1U|) , for single-inclusive neutral pion production which 
will be one of the main quantities of interest in experiment. 
Fig. E| shows A^f, calculated at NLO for the "standard" 
set of the GRSV spin-dependent parton densities |2"4"|. as 
well as for three other sets emerging from the GRSV anal- 
ysis which mainly differ in the assumptions about Ag (see 
above). The impact of the unknown non-perturbative par- 
ton structure of the circularly polarized photon on A^jf is 
examined by making use of the two extreme sets also intro- 
duced at the beginning of Sec. 01 We refrain from showing 
LO estimates for the double-spin asymmetry which are of 
rather limited use anyway. Due to the pronounced differ- 
ences in the if-factors for the unpolarized and polarized 
cross sections, see Fig.|2 the LO results for A££ are con- 
siderably larger than the NLO ones shown in Fig.0J This is 
in contrast to the frequently made assumption that NLO 
corrections cancel in spin asymmetries. We note that simi- 
lar observations have been also made for hadroproduction 
of pions [S]. 

As can be seen, the actual choice of photonic parton 
densities barely affects the results for the spin asymme- 
try shown in Fig. 0] if the pion's transverse momentum is 
larger than about 2 GeV. This can be readily understood 
by noticing that in this region the average momentum 
fraction (x a ) in Eq. J5J is larger than 0.5, i.e., one probes 
only x 7 -values where the photon structure is dominated 
by the "pointlike" contribution independent of the details 
of the unknown non-perturbative input |28j . In this px- 



region a measurement of the spin asymmetry could be re- 
lated to a certain Ag. However, this does not imply that 
the resolved contribution to the cross section is negligible 
as we shall demonstrate below. On the other hand, for 
Pt ;$ 1.5 GeV results for different gluon polarizations of 
the nucleon and assumptions about Ap overlap, making 
it virtually impossible to draw any conclusions without 
either knowing Ag or Af 1 . 

To judge whether the observed dependence of A^jf on 
Ag in Fig. g| can be used to learn something about the 
gluon polarization, we give estimates for the expected sta- 
tistical accuracy SA^ for such a measurement at COM- 
PASS in certain bins of px, calculated from 



5A 



H,N 
LL 



1 



1 



(11) 



Here, Ubin denotes the unpolarized cross section in the px- 
bin considered, C the integrated luminosity for which we 
assume 1 fb~ , and all other parameters are as specified at 
the beginning of Sec. 13.11 Not unexpectedly, the statisti- 
cal accuracy rapidly deteriorates towards higher px as the 
photoproduction cross section drops sharply. In the region 
1-5 < Pt S 2.5 GeV it is conceivable that some informa- 
tion about Ag can be obtained. We note that the gluon 
polarization will be predominantly probed in the range 
0.1 < .Tf, < 0.3 of the momentum fraction Xb in Eq. (J2J. 

Let us now turn to a closer analysis of how the re- 
sults presented in Figs. |3 and 0] come about and can be 
understood, by studying the different contributions to the 
polarized photoproduction cross section separately. This is 
done in Fig.JSJ Here we use again our default sets of parton 
densities: the "standard" set of GRSV [21] for the nucleon 
and the "maximal" set of for the photon. First of all, 
Fig. El (a) reveals that the resolved photon cross section 
dAa rcs dominates over the direct contribution <iZ\<7di r in 
the entire range of px relevant for Compass. At a first 
glance this result is counter-intuitive and requires further 
clarification since at fixed-target energies one might ex- 
pect a dominance of the direct contributions to the full 
photoproduction cross section. The importance of the re- 
solved photon contribution, in particular in the polarized 
case, can be readily understood by a closer inspection of 
the individual partonic subprocesses contributing to the 
full photoproduction cross section. 

Figure|5](b) shows that indeed both direct channels, ■yq 
and jg scattering, exceed in magnitude any of the resolved 
subprocesses. However, the two channels contribute with 
different sign and therefore cancel each other to a large 
extent in the sum dAa^h-. In this way numerically smaller 
resolved subprocesses can become important, in particu- 
lar, as they all contribute with the same sign for a positive 
gluon polarization. We note that the partonic spin asym- 
metry dAa/da for the LO process -fg — ► qq is —1, such 
that only positive Ag lead to this cancellation pattern. 
Nevertheless, also for a large negative gluon polarization 
["Ag = —g input" in Fig. 0(b)] is the resolved contribu- 
tion non-negligible. For all relevant values of px the bulk 
of the resolved contribution stems from the qg subprocess 
where a quark at large momentum fraction x 7 from the 
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Fig. 5. (a) Direct and resolved NLO(MS) contributions and 
(b) NLO(MS) contributions of different partonic channels a + 
b — > c + X to the full NLO polarized photoproduction cross 
section shown in Fig. 



photon scatters off a gluon from the nucleon [solid line in 
Fig. El (b)]. This is because quark densities in the photon 
are sizable at large x 1 |2.'iU28| due to the pointlike contri- 
bution, contrary to parton densities in the nucleon which 
rapidly vanish as x — » 1. As mentioned before, at large 
momentum fractions Aq 1 is not sensitive to the details of 
how the non-perturbative hadronic content of the photon 
was modeled. We therefore conclude that photoproduc- 
tion of inclusive hadrons with Compass kinematics and 
statistics will most likely not yield much information on 
the non-perturbative nature of the resolved photon. On 
the other hand, this finding certainly simplifies an analy- 
sis of photoproduction data with px > 1.5 GeV in terms of 
Ag despite our ignorance of the photonic parton densities 
Ap. 

We obtain very similar results when considering a pro- 
ton target, as will be realized when the COMPASS experi- 
ment switches to a NH3 target in the future. The resulting 
spin asymmetry, A^, is depicted in Fig. To estimate 
the statistical accuracy we have again assumed an inte- 
grated luminosity of lfb -1 , and all other parameters are 
as specified at the beginning of Sec. 13.11 The sensitivity 
to Ag and the rather weak dependence on the photon 
scenario for px > 1.5 GeV, characteristic for the large x 7 - 
region probed at COMPASS, are essentially the same as for 



■ We note that the intricate interplay between direct 
and resolved photon processes is quite similar to the one 
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Fig. 6. As in Fig. [I] but now for a "proton" target. 
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Fig. 7. As in Figs. 2] and 03 but now for the photoproduction 
of charged hadrons (see text). 



described in Figs. 0(a) and (b), but that the cancellation 
between the 73 and 75 channels is less complete such that 
the resolved cross section is somewhat less relevant here. 

So far we have only considered the production of neu- 
tral pions where fragmentation functions were shown to 
work reasonably well also at rather low scales How- 
ever, the sum of charged hadrons, predominantly pions, 
but also kaons and protons, is equally important exper- 
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Fig. 8. Unpolarized and polarized p-r- differential cross sections 
at NLO for the reaction /id — > ^l-k°X for two different exper- 
imental setups: 9 mal = 180 mrad (solid) and, as in Fig. [5] 
flmax = 70 mrad (dashed). The lower panel shows the corre- 
sponding ratios of NLO and LO results (AT- factor). 
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Fig. 9. As in Fig. 0] but now for max = 180 mrad. 



imentally as these are often more easily identified than 
neutral pions. This is also the case for Compass at the 
moment. In Fig. we therefore present the relevant spin 
asymmetries A^l and A^f f° r deuteron and proton tar- 
gets, respectively, for the reaction fid(p) — > fi'hX, where 
h represents the sum of charged hadrons (pions, kaons, 
and protons). The results are obtained by employing the 



appropriate set of fragmentation functions, +fl of 
Rcf. [23. Not unexpectedly, the gross features of the spin 
asymmetries in charged hadron production are the same 
as in neutral pion production. Due to the larger rate for 
the sum of charged hadrons, the statistical precision is 
noticeably better than for A£jf which, makes such a mea- 
surement potentially more interesting. 

From our results for the spin asymmetries shown in 
Figs. El and it is obvious that a major difficulty in 
extracting the gluon polarization Ag at fixed-target ener- 
gies lies in the poor statistical accuracy at large p^-values. 
This, however, is the region where perturbation theory is 
expected to be more reliable and where the uncertainties 
associated with the resolved photon contributions to the 
cross section are much better under control. It is there- 
fore interesting to study whether the planned upgrade of 
the Compass experimental setup, which would lead to a 
much larger acceptance of max = 180 mrad and hence 
larger cross sections, could help. 

In Fig. |H1 we compare the polarized and unpolarized 
neutral pion photoproduction cross sections for the present 
setup with # max = 70 mrad and for the proposed upgrade 
with # max = 180 mrad. The gain in cross section is a factor 
4-^5 depending on the pr-value considered, yielding an 
improvement of the statistical accuracies for spin asym- 
metry measurements by about a factor of two. The lower 
panel of Fig. |S| shows that the respective if- factors do not 
change significantly when going from # max = 70 mrad to 
#max = 180 mrad. However, not only does an increased an- 
gular acceptance improve the statistical accuracy, it also 
changes the interplay of direct and resolved subprocesses. 
It turns out that differences between the "maximal" and 
"minimal" Af 1 scenarios are more pronounced and persist 
towards larger px- This can be inferred from our estimate 
for the spin asymmetry A^jf shown in Fig. if compared 
to the results given in Fig. Also the sensitivity of A£ \ 
to different assumptions about the gluon polarization of 
the nucleon is somewhat smaller than before. 



3.2 Single-Inclusive Hadron Production at HERMES 

Single-inclusive hadron photoproduction can be also stud- 
ied with the Hermes experiment at DESY where an elec- 
tron (or positron) beam with E e ~ 27.5 GeV is scattered 
off a proton or deuterium gas target. The available c.m.s. 
energy of about — 7.5 GeV is lower than at COMPASS 
which even further limits the range of accessible trans- 
verse momenta. On average the electron beam polariza- 
tion is V e ~ 53%. For the polarization of the gas target 
we take Vd ~ V v — 85%, and, contrary to a solid-state 
target, there is no dilution, i.e., T p = Td = 1. To esti- 
mate the statistical accuracies for a measurement of A^^ 
and A£l* below, we assume an integrated luminosity C of 
50 pb -1 and 200 pb -1 for proton and deuterium targets, 
respectively. These numbers are based on the actual data 
already collected with the Hermes experiment which is 
now running with transverse polarization most likely un- 
til the end of their experimental program. 
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photoproduction of neutral pions, ep 
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The lower panel shows the ratios of NLO to LO contributions 
(if- factor). 
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Fig. 11. As in Fig.[3](b) but now for ep — » e'-K°X at Hermes. 



The acceptance of the Hermes experiment allows the 
detection of pions in the angular range 40 mrad < 9 < 
220 mrad. The upper limit on 9 corresponds to a mini- 
mal c.m.s. rapidity of r;™" — 0.2. At small pr the lower 
limit on 0, i.e., i]™* — 1.91, sets a more stringent bound 
on r]™* than purely kinematical considerations, f?™^* = 
cosh" 1 ( % /5/2p T ). For all our numerical studies we inte- 
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Fig. 12. Double-spin asymmetries A^jf and A^' L f° r Hermes 
kinematics at NLO for different gluon polarizations in the nu- 
cleon and minimal and maximal saturation of the polarized 
photon densities, dashed and solid lines, respectively. The "er- 
ror bars" indicate the estimated statistical uncertainty assum- 
ing an integrated luminosity of 50 pb -1 and 200 pb -1 for -A^jf 
and > respectively. 



grate over rapidity as before. We select photoproduction 
by demanding that the virtuality squared of the photon is 
less than <5„ lax = 0.01 GeV 2 . The photon flux with pho- 
ton energies limited to the range 0.2 < y < 0.9 is again 
modeled by Eq. @. 

The resulting polarized and unpolarized px-diffcrcntial 
cross sections and AT-factors for photoproduction of neu- 
tral pions are depicted in Fig.^J We only show results for 
a Hermes setup with a proton target, i.e., ep — > e'w°X, 
as results for a deuterium target are qualitatively very 
similar. Due to the small c.m.s. energy available the ex- 
perimentally accessible transverse momenta pt are limited 
to a region where a successful and reliable application of 
perturbative methods cannot be taken for granted. This 
is emphasized, for instance, by the large value of the un- 
polarized A"-factor which ranges from 2 to 4. Even more 
dramatic are the uncertainties related to the choice for the 
renormalization and factorization scales which are shown 
for the polarized photoproduction cross section in Fig. ITT1 
Again we have varied all scales simultaneously in the range 
Pt/2 < fJt r = fj,f = fjft < 2pr- There is only a slight im- 
provement when the NLO QCD corrections are taken into 
account. 
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Fig. 13. (a) Direct and resolved NLO(MS) and (b) NLO(MS) 
contributions of different partonic channels a + b ^ c + X to 
the full NLO polarized photoproduction cross section shown in 

Fig. ina 



These results again clearly indicate that special care 
has to be taken in the interpretation of upcoming results 
in terms of pQCD and, ultimately, of Ag at current fixed- 
target experiments. It is perhaps even more important 
than at Compass energies that studies of the theoret- 
ically better predicted unpolarized cross section accom- 
pany studies with polarization. All reservations notwith- 
standing, we analyze in Fig. 1121 the sensitivity of hadron 
photoproduction at Hermes to the gluon polarization of 
the nucleoli. We show the double spin asymmetries 

and a ^ NLO f° r different assumptions about the gluon 
polarization of the nucleon and the polarized parton densi- 
ties of the photon. Again, we estimate the statistical error 
for such a measurement with the help of Eq. I|ll|) with 
the parameters as specified at the beginning of this sec- 
tion. Clearly, meaningful results about Ag can be only 
obtained for px ~ 1.5 GeV. Due to the larger data sample 
available, results obtained with the deuteron target are 
much more discriminating than with the proton target. 

Figure 1721 also shows that the uncertainties from esti- 
mating the resolved contribution to the cross section are 
not very pronounced even at the lowest j>T-values shown. 
To examine the importance of the resolved contribution at 
Hermes further, we show in Figs. El ( a ) an d (b) the rela- 
tive contributions of the direct and resolved and different 
partonic subprocesses a + b —>■ c + X to the full photopro- 



duction cross section, respectively. It turns out that for a 
proton target the direct photon contribution dAassx dom- 
inates in the entire p^-range. Contrary to the situation for 
COMPASS, cf. Figs. GO (a) and (b), the cancellation of the 
two direct channels 73 and 79 in the sum dAadir is less 
complete. All resolved subprocesses are small and, for the 
same reasons as discussed before for Compass, the bulk 
of dAa rcs stems from the scattering of a quark carrying a 
large momentum fraction of the photon off a gluon in the 
nucleon [solid line in Fig. ^5] (b)]. For a deuteron target 
one obtains qualitatively similar results. However, at the 
small pr-values shown the resolved contribution is more 
important due to a more pronounced cancellation of the 
75 and jq subprocesses in dAadir- 

3.3 Comparison to the E155 Data 

As we mentioned in the introduction, there is already a 
data set on the spin asymmetry in single-inclusive hadron 
production by the El 55 collaboration |1G[. and it is of 
course interesting to compare our NLO calculations to 
these data. The main data set of E155 is for eN — > H ± X, 
where H + (H~) denotes any hadron of identified positive 
(negative) charge. Only for the data at the larger scat- 
tering angle 9 — 5.5° are the observed hadron transverse 
momenta large enough for a sensible comparison to pQCD 
hard-scattering. Here we use again the various GRSV sets 
of polarized parton distributions, and the "maximal" set 
of polarized photon densities for the resolved contribution. 
We use the fragmentation functions of Ref. [3U| which pro- 
vide separate sets of H + and H~ fragmentation functions. 
Finally, because the scattered electron was not observed 
for the El 55 data, we use the Weizsacker- Williams spec- 
trum in Eq. © in the case when all photon virtualities 
are integrated, and we do not impose a cut on the variable 

Figure ITU shows the comparison of our NLO calcula- 
tions to the data of E155 for H + and H~ production 
off proton and deuteron targets. One can see that it is 
difficult to achieve a quantitative overall agreement for 
a given set of GRSV polarized parton densities. On the 
other hand, the results show very large variation with the 
polarized parton distributions, and fine details of the den- 
sities will therefore matter. This is in particular so as there 
are again strong cancellations between the direct and re- 
solved contributions, and even within these. One should 
also keep in mind that the uncertainty in the charge sep- 
aration of the available sets of H + and H~ fragmenta- 
tion functions is currently very hard to quantify and cer- 
tainly large. We finally note that there is also room for 
the "soft" non-perturbative contributions to the asymme- 
try discussed in even though the data do not give 
compelling evidence for their presence either. 

4 Discussion and Conclusions 

To summarize, in this paper we have presented a detailed 
phenomenological study of single-inclusive hadron photo- 
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Fig. 14. Comparison of our NLO calculations to the E155 
data for the double-spin asymmetries for H + and H~ produc- 
tion in electron scattering off proton and deuteron targets at 
scattering angle 9 = 5.5° |16l . In each of the four panels, differ- 
ent line styles distinguish the different GRSV sets of polarized 
parton distributions used in the calculation. 



production at fixed-target energies with particular em- 
phasis on the kinematics relevant for the COMPASS and 
Hermes experiments. We have also compared our calcu- 
lations to already available data from the El 55 experi- 
ment. All calculations were consistently performed at the 
next-to-leading order of perturbative QCD. We have care- 
fully explored and critically discussed the applicability of 
the perturbative methods in the region of relatively low 
momentum transfers characterized by the pT-values ac- 
cessible in fixed-target experiments. Both the if-factors 
as well as uncertainties due to the choice of the renormal- 
ization and factorization scales turn out to be sizable; the 
latter only slightly improving when NLO corrections are 
taken into account. 

All this indicates that the pQCD framework in the 
Compass and Hermes kinematic regimes is delicate, and 
calls for a "reference measurement" of the unpolarized 
photoproduction cross section before the spin asymmetry 
can be analyzed in terms of Ag. If the unpolarized cross 
sections turn out to be in reasonable agreement with the- 
ory, measurements of the corresponding spin asymmetries 
by Compass and Hermes will add valuable information 
about the gluon polarization of the nucleon at momentum 
fractions x ~ 0.1 0.3. To access these large a;- values at 
a polarized pp collider like RHIC requires to detect final- 
state hadrons, photons, or jets at rather high px where the 
rate is small. The theoretical analysis of the unpolarized 
cross sections at Compass and Hermes may well require 
to pursue avenues like performing all-order resummations 
of large logarithms in the perturbative series. 



Finally, it turned out that the resolved photon contri- 
bution to the cross section is sizable, and, under certain 
conditions, even dominant, but independent of the details 
of the non-perturbative spin structure of the polarized 
photon which is completely unknown so far. Analyzing 
the data with the approximation dAa ICS would lead 
to incorrect conclusions about Ag. To access the parton 
content of polarized photons certainly requires a future po- 
larized lepton-hadron collider facility such as the eRHIC 
project at BNL. 
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